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Abstract - In this paper, we consider a model of interacting phantom and Yang-MiUs 
(YM) fields supposing the dilaton-type coupling. Making use of specific solution for 
YM equation previously found by the author, we obtain simple exact solutions for the 
, accelerated expansion of FRW cosmological model. Besides, we derive the induced 

04 ' potentials of phantom field corresponding to some given regimes of expansion. The ef- 

fective equation of state (EoS) have been reconstructed for all types of the models 
^ , considered here. 
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O^' 1 Introduction 

u . 

Present accelerated expansion of the universe is well proved in many papers [l]-[8]. In order 
to explain so unexpected behavior of our universe, one can modify the gravitational theory 
CO ' or construct various field models of so-called dark energy (DE) which EoS satisfies 

w = p/ p < —1/3. The simplest candidate of DE is the cosmological constant with fixed 
EOS w = —I. If it is quintessence then — 1 < ui < — 1/3 and if it is phantom then w < —1. 
The constant EOS w — —1 is called phantom divide. There are some dark energies which 
^SJ ' can cross the phantom divide from both sides [Ij]. So far, a large class of scalar- field DE 

models have been studied, including tachyon [TB], ghost condensate [T7] and quintom [TB], 
' [12], and so forth. In addition, other proposals on DE include interacting DE models |20| . 

braneworld models "IV , and holographic DE models [35], etc. The quintom scenario of DE 
is designed to understand the nature of DE with EoS across -1. The quintom models of DE 
differ from the quintessence, phantom and k-essence and so on in the determination of the 
cosmological evolution. It is appropriate mention here that impossibility to realize quintom 
(or crossing of the phantom divide) in genera k-essence-like models was shown in [23^. 
' Another class of DE models is based on the conjecture that a vector field can be the 

origin of DE [24], [25]. The YM field can be a kind of candidate for such a vector field [26], 
[27] . [28) . At the same time, it is well known that a pure YM field (with its EoS w = 1/3) 
can not provide accelerated expansion of the Universe, for which w < —1/3 is required. 
This is a direct consequence of conformal symmetry of the Lagrangian for a massless YM 
field. Any violations of conformal symmetry (e.g., as a result of quantum corrections [29] 
or of non-minimal coupling to gravity 3D]) give a good chance for involving YM fields in 
reconstruction of DE. The alternative chance for YM fields to be involved in DE problem is 
consideration of some interaction of YM field with different sources of gravity. In this aspect, 
the idea of induced nonlinearity is fairly attractive for realization in cosmology of phantom 
field. In this paper, we turn our attention to the issue of the YM fields interacting with a 
phantom field in FRW cosmology. Using the specific solution of YM equation previously 
considered in FRW cosmology [31]- [35], we generalize the model investigated in [31] on the 
case of interacting phantom and YM fields. This allows us to obtain some exact solutions 
for the accelerated expansion of FRW cosmological model. Besides, we derive the induced 
potentials of phantom field corresponding to some given regimes of expansion. The effective 
EoS have been reconstructed for all types of the model considered below. 
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2 Basic equations 

The main equations of the model follow from the Lagrangian density [31j : 

+ (1) 

where R is the Ricci curvature scalar, 1^9 is a scalar field, = daW^ — dpW^ + geabcW^i^Wp 
is the YM strength tensor, ^ (<p) is the coupling analytical function of the phantom and 
YM fields, e = +1 represents quintessence while e = — 1 refers to phantom field. 

As it is noted in [31], the similar Lagrangian is well founded by the existence of pro- 
cess TT — >■ 27. The case ^((p) = exp(— 2A(/?) with A — \/3 occurs via the Kaluza- Klein 
compactification of five-dimensional vacuum gravity. The critical coupling A = 1 arises by 
the truncation of = 4 supergravity.The same type of interaction emerges in Brans-Dicke 
theory with A — \/\/2u] + 3. This fact encourages us to suppose that the theory ([T]) can be 
embedded into some cosmological theories for general coupling. 

Variation of dD) with respect to metrics g'^'' and fields yields the Einstein equation 

G", - ^t;, (2) 

where GJ^ is the Einstein tensor, k is the Einstein gravitational constant. From ([1]), one can 
obtain the following energy-momentum tensor for the system of fields, 



and the modified YM equation: 



(3) 



{^gF'^-^^{^)) = 0, (4) 
where Di, denotes the covariant derivative. The scalar field equation is as follows: 



5.9^^7^ +71^^*. = 0, (5) 



—g dx" \ ' dx>^ J IGtt 



where = — and the first invariant of Yang-Mills field / = F^^F'"^'^. 



dip 

We assume that the Universe is described by a Friedmann-Robertson- Walker (FRW) 
geometry: 

ds^ = df - a^{t){dr^ + i^{r)dn^), (6) 

where ^(r) = sinr, r, sinhr for the sign of space curvature k = -1-1,0,— 1, consequently. As 
well-known, the generalized Wu-Yang ansatz for the YM fields can be written as [3T] 

gr ^ gr'^ V r / gr 

We can make the following substitution into this ansatz |33j : 

W{r,t)=a{t), K{r,t) ^ P{r) cos a{t), S{r,t) = P{r) sma{t). 
As a result, we have the following expressions for the YM strength tensor components: 

Foi = F02 = Fo3 = 0, F12 = g^^P'ir)(in cosa -f 1 sina^, 
F13 = g^^P'{r) sin sin a — 1 cos ct^ , 

F23=5"isin0(p2(r)-l)n, (7) 

where n = (sin 61 cos 0, sin sin 0, cos 6*), 1 = (cos 6* cos 0, cos 6* sin 0, — sin 0) and m — 
(— sin 0, cos 0, 0) are the orthonormalized isoframe vectors, and the prime means a derivative 
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with respect to r. As noted in [33], the YM field (O has only magnetic components. It is 
easy to find from ^ and (0 that the YM field invariant / = F°-f,F'^'^^ becomes as follows: 



/ = 



(8) 



With the help of ([7]), assuming the spacial homogeneity of phantom field , it is easy to show 
that YM equation Q reduces to the following simple equation |31) : 

fp2 _ DP 

P" ^ = 0' (9) 

where P" = (fP/dr^. 

The nontrivial solution for equation © obtained in 31 is as follows: P{r) — — 
COST, cosh r for k = +1,-1, consequently. From the latter and ([5]), it follows that the 
valuable feature of this solution is that the YM invariant built on this solution depends only 
on time: 

^ = ^W = ^477T- (10) 
In view of ([TU|) . the rest set of Einstein and scalar field equations is as follows: 

-Ja' + k) = ne^ + -^^n^), (11) 



'"^ 2aa + fc) = ~ 
15,,.. 3 1 



2(a2 + 2aa + fc) = -«e^--^^*(^), (12) 



It is easy to verify that only two equations among Eqs. (|11I) - P^ are independent. That 
is why we should define one of the three functions a{t), ip{t) or ^'(95) to solve the system. 
The most natural way is to define the coupling function. But it is not necessary and depends 
on a specific context of the problem. In works |31| , |33| , this set of equations is investigated 
from the assumption of some given regimes for expansion, fc = ±1 and e +1. Now we 
consider the phantom case, that is e = —1. 



3 Simple examples of exact solution 

Let us study the set of dynamics equations ([TT |) - (|13p for the case of phantom field consid- 
ering e = —1. The set of Eas. pT|) - p^ can be re-written as the following two independent 
equations for ip and ^'(((s): 

^'^ ^Gn^'^a^^ia'^ + ad + k), (ip) ^ Sirg'^K^^a^ {2d'^ + ad + 2k). (14) 



Moreover, comparing Eq. (jl3p in the case e — —1 with the usual one for a phantom field in 
FRW cosmology, 

dip 

where H = — is the Hubble parameter, we derive the relationship between an effective 
a 

potential V^ff and the coupling function ^'(1^9): 

dV^ffjifi) _ 3 1 d^jip) ^^^^ 
dip 87rf/2 a4 

From the first equation of (|14p it follows that the real solution for phantom field exists 
when 

+ ad + k>0. (16) 
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Besides, the simple study of effective energy density and pressure which foUow from the 
r.h.s. of (dll) and ([n]), 



87rg2a4 



(17) 



shows that the effective EoS in this model is restricted by p < — p/3. Therefor, the acceler- 
ated expansion, for which p < — p/3 is required, can be realized in our model. 

To demonstrate some interesting features of this model, we are going to consider two 
illustrative examples of exact solutions further to our general study. We will consider the 
accelerated expansion, that is a > 0, when the asymptotic value of the Hubble parameter 
H = Hq — const is achieved in the course of time. For all that, the necessary condition (ITBl) 
is satisfied at every instant. 



3.1 The case a{t) = H^^ smh{Hot) 

Note that this dependence of a{t) on time is satisfies inequality ([T5)) for all signs of the 
curvature, and Eq. (fT4)) becomes as follows: 



6HS 



K sinh^ (Hat) 
Sng'^ sinh^(iJot) 



(2smh^{Hot) + l + k 



(3 sinh^ (Hot) + 2 + 2k 



(18) 



(19) 



Consider the cases of negative and positive signs of curvature separately. Here and further, 
we use the following notations: 



A 



24Trg^ 



Vo 



3Hl 
A2 



i) Open model: A; = — 1. From Eqs. ([T5| . (IT^ . one can find that 

ip = X-^Ho, l'(t) = 5-0 sinh'*(i/oi)- 



(20) 



Integrating the first equation and taking into account the second one together with the 
explicit expression of a(<), we obtain that: 



ip{t) = X ^Hot + ipo, = ^osinh"* X{ip - ipo) 



(21) 



With the help of equation dTSl) . we can derive the following effective potential which is 
induced by the coupling to the YM field: 



Veff{ip) = Vo In (^sinh X{(p ~ ipa) ^ + Uq, 



(22) 



where Uo is a constant. 

ii) Closed model: k = +1. As it follows from Eqs. (IT51) . ([T^ in this case, we have 



<p = X-^HoCoth{Hot), *(<) = ■i>osmh'^{Hot)(smh^{Hot) + 
Then from Eq. 



(23) 



we can find that 



ip{t) ~ X sinh(_ffoi) +^0, 



(24) 



The induced phantom potential, which corresponds to the coupling function (|24p . can be 
obtained in the following form: 



Veffiv) = Vo X{(p - (po) 



1 -2A((^ - </.o) 



+ Uo. 



(25) 
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Figure 1: Effective potential 14// {Vq — 1,Uo ~ 0) versus f — for a ~ H^^^ amh{Hot) 
when fc = — 1 (dashed curve) or k = +1 (dotted curve), and for a = Hq^ cosh(iJoi) when 
k = —1 (chain curve) or fc = +1 (sohd curve). 

3.2 The case a{t) = H^^ cosh{Hot) 
In this case, Eqs. (11^ take the foUowing form: 



Kcosh (-ffo^) 



(2cosh^{Hot)-l + k 



87r.g2 cosh^(i?oi) 



(^3cosh^(i7ot)-2 + 2fc), 



iii) Open model: fc = — 1. From Eqs. (pS)) . (P7)) we have 

= A"^iJotanh(iJot), *(i) = *o cosh^(i?oO(cosh^(iJoO - - 

V 3 



(26) 
(27) 

(28) 



In view of the second equation in (j28p and explicit expression for a(t), we can solve the first 
equation in (|28p . and then obtain 



Lp{t) = A"^ In cosh(i/oi) 



*(^) = ^'oe2^(^ - ^o) [e2'^(^ - ^o) - - 

3 ■ 

Due to Eq. ([T5|) . the effective phantom potential becomes 

f 



(29) 



y,//(^) = [a(^ - ^o) + -e-2A(^ - (^o) 



C/o- (30) 

iv) Closed model: fc = +1. For this sign of curvature, it is follows from Eqs. ([2^ . ([57]) that 

ip = X-^HQ, l'(t) = ^'ocosh*(i7oi). (31) 
Therefor, we have the following solution for (|3ip: 

(^(i) = A-iiJoi + (^o, ^'(^) = 1-0 cosh^ [a((^ - v?o)] • (32) 
In view of (|15p . it can be easily obtained that now the effective phantom potential is 

K//((^) = 14)ln(cosh[A((^-¥>o)]) +Uo. (33) 



The plots of effective potentials for all cases considered can be viewed on Fig.l. 




Figure 2: Effective EoS parameter w versus x = Hot for a = H^^^ sinh{HQt) when k = — 1 
(dotted curve) or fc = +1 (solid curve), and for a = Hq^ cosh(i?oi) when fc = — 1 (dashed 
curve) or fc = +1 (dotted curve). 

With the help of expressions ([T7|. all solutions obtained above can be arranged in tree 
groups according to their EoS. 

In case (1): a(t) = Hq^ cosh{Hot), (fc = —1), and in case (2): a{t) = Hq^ smh{Hat), 
(fc = +1), one can find consequently 



Wi 



3smh^{Hot) + 1 
3sinh2(i?oi) -3' 



W2 



3sinh^(iJoi) +2 
3sinh2(iJot) +6 



(34) 



In case (3): a{t) ~ H^^ sinh(iJoi) with fc —1, or a{i) ~ Hq^ cosh{Hot) with fc — +1, we 
have ^3 = —1. The plots of EoS parameters for all three cases considered above can be 
viewed on Fig. 2. 

As one can see, the curve of EoS in the case a — Hq^ sinh(i/oi) with fc = — 1 displays 
rather unusual behavior around the value x = Hgt determined by the solution of sinh(iJot) = 
1. To understand this, let us find the energy densities from Eq. ([T7|) for the cases defined 
above. So we can find for the cases (1) and (2) that 



Pi 



4A2 cosh''(iJo<) ^ 



( sinh^(iJot) - 1 



P2 



Hi 



4A2 smh^H^t) V 



( sinh^(iJot) + 2 



(35) 



(36) 



For case (3), we have ps 



Hi 



, that is p3 = lim pi (t) = lim p2 (t) . It could be supposed 

4A t— >00 t— >CXD 

that the constant density pa represents some energy density of an effective cosmological 
A 

constant A, that is — = — ^r. This implies the well-known result for a de Sitter model: 
K 4A2 

Hq — K/i. At the same time, one can conclude from ((35)) that while smYv' {H^t) < 1 the 
energy density pi is negative, i.e. the weak energy condition is violated. But as soon as 
sm^c? {HqI) > 1, pi becomes positive, and the energy condition takes effect. From (|36|) . it is 
obvious that there is no any problem of this sort in case (2), as well as in case (3). 



4 Conclusion 

In summary, the model of interacting phantom and YM fields in FRW non-fiat cosmology 
are shortly studied in this paper. First of all, we have derived the set of main equations 
which determines the model dynamics: (fTTj) . ([T^ and (fT5|) . Making use of the specific 
solution of YM equation previously considered in FRW scalar field cosmology, we generalized 
the model investigated in [29] on the case of interacting phantom and YM fields. This 
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allowed us to obtain some exact solutions for the accelerated expansion of FRW cosmological 
model. Besides, we derive the induced potentials of phantom field corresponding to the 
cases (3.1) and (3.2) in which the Hubble parameter changes as H{t) = iJp coth(i7oi) or 
H{t) = i?o tanh(i?oi) consequently. At that, it is follows from ^ and ^ that 

all cases considered above are related by the conditions </5 ^ Hq or (ys ~ H{t). In other 
words, the rate of phantom field change is proportional to either asymptotical value of the 
Hubble parameter ( Hq — lim H{t) ) or to its contemporary value. The effective EoS have 

been reconstructed for all types of the models considered above. Somewhat unexpected is the 
results of cases (1). Nevertheless, this model considered from the moment t — Hq^ sinh~"'^(l) 
does not demonstrate any oddity. As can be seen from our examples, all EoS parameters 
considered do not cross the phantom divide —1 at late time. This is not a characteristic 
property of the model but only the consequence of the simplest expansion regimes given. 
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